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Abstract 

The equations of viscous evolution of 3D arbitrarily shaped vortices in 
an isotropic type II superconductor and necessary boundary conditions are 
formulated in the frame of London approximation. The theory is applied to 
analyse characteristic scenaria of vortex penetration into current-carrying 
thick plate or bulk samples with another geometry. It is shown that re- 
garding of the surface transport current value a vortex penetrates either as 
"flexible stick" or as "elastic thread". The latter regime is accompanied by 
giant stretching of the vortex core along the current-induced surface mag- 
netic field. This geometrical transformation leads to decrease of viscous 
friction and large increase of the vortex drift velocity as compared with 
the stick-like regime. As a result, the vortex first winds round the sample 
cross-section and forms a ring-like curve, and only later begins to move deep 
into the sample interior. The analytical estimates of the vortex shape and 
stretching and its velocity are obtained. 

I. INTRODUCTION 

Not long ago a variety of complex magnetic structures formed by many strongly curved 
and entangled vortices was discovered in bulk superconductors [1]. The origin of these 
structures can not be explained if treat the motion of vortices like that of stick-like objects. 
It is necessary to consider the evolution of three-dimensional magnetic flux lines with 
potentially arbitrary shape. 

For the first, it is useful to investigate the motion and shaping of a single vortex 
but interacting with a surface supercurrent which represents either transport current or 
Meissner current induced by external field. As far as we know, even this simple task 
previously was not under careful consideration. 

Of course, a vortex never lives alone, without interactions with other vortices, and no 
stable many-vortex structure could exist without mutual repulsion of vortices. However, 
one can suppose that the scenario of magnetic flux penetration into a current-carrying 
bulk superconductor should be dominated only by vortex interaction with the surface 
current distribution, i.e. eventually by geometry of superconducting sample, not by the 
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vortex-vortex interaction. The latter can not seriously affect this scenario, merely because 
it itself is unable to ensure a deep penetration at all. 

To prove this statement, let us imagine the steady flow of vortices which arise at the 
flat current-carrying boundary of half-infinite superconductor and then move deep into 
the sample due to their repulsion. Clearly, because of viscous character of vortex motion, 
such the flow needs in nonzero gradient of concentration of vortices. As a consequence, 
both the concentration and the local drift velocity of vortices must be decreasing functions 
of the depth. Hence, their product is not constant, that is the magnetic flux conservation 
can not be satisfied. This discrepancy means that no steady flow could be supported by 
the inter- vortex forces only. In particular, it is impossible to realize the stationary lasting 
penetration of vortices from infinite flat boundary parallel to external magnetic field. 

Therefore, the only force what can push a vortex through the sample interior is nothing 
but self-action of vortex caused by its distorsion. But in order to involve this force into 
the evolution, the vortex must feel the shape of the sample boundary. 

Hence, the true picture looks as follows. After nucleation in a surface layer with 
thickness of order of London penetration depth A , a vortex firstly expands over the 
sample boundary remaining in this layer. At this stage only the end fragments of vortex 
are factually moving. The ends slide along the boundary, and the resulting shape of vortex 
core reflects that of boundary. This process lasts until the curvature of main middle part 
of the core becomes sufficiently strong in order to cause the deepening of vortex as a 
whole. 

In view of these reasonings, the geometry of steady transport of magnetic flux into a 
bulk supeconductor looks rather insensible to inter-vortex interactions and thus can be 
testified in terms of a single vortex, at least if not consider details of vortex nucleation 
and processes like annihilation and reconnections of vortices [1] which take place deep 
inside the sample. 

For example, many aspects of resistivity in supercurrent-carrying wires can be de- 
scribed as evolution of ring-like vortices as if thats instantly arise near the boundary, 
then contract independently one on another and finally self-annihilate [2]. However, more 
correct consideration should include the first stage when vortex transforms from small 
nucleus into a ring. We shall see that in fact this stage may result also in a non-ring 
penetration geometry, and more detailed theory can predict what the scenario realizes 
under given transport current value and sample dimensions. 

Though a lot of works were published previously touching upon a role of vortex dis- 
torsions, for instance, under a pinning by randomly distributed centers [3], always some 
preliminary restrictions of the vortex geometry were attracted. In the present work the 
general equations of evolution of arbitrarily curved vortex lines in isotropic superconduc- 
tors are formulated and analysed. We shall especially discuss the true formulation of 
boundary conditions for these equations. 

It will be shown that in a sample whose dimensions noticably exceed A the vortex can 
penetrate either as flexible stick or as elastic (similarly to a rubber thread attracted by 
its ends through water), of surface supercurrent. The latter case occurs only if surface 
supercurrcnt exceeds Hcic/4tt (in CGS units) and is characterized by giant stretching 
of the vortex core along the sample boundary in the direction parallel to drift of the 
ends. The stretching is accompanied by decrease of both the vortex energy and viscous 
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dissipation per unit drift velocity, and results in strong increase of the vortex drift velocity 
under given transport current. In the framework of this scenario, the vortex core firstly 
tranforms into a ring-like curve winding round the wire cross-section (or into a spiral, if 
there is an external magnetic field parallel to current), and only later the vortex begins 
to cut the wire and enter deep into its interior. This general picture is in agreement with 
the known simplific model of magnetic fiux penetration into round wires. Additionally, 
our approach allows to scope very different stages of vortex evolution in unified manner 
and obtain quantitative estimates for each stage. 



II. LONDON APPROXIMATION 



We shall confine ourselves by the London approximation. Of course, it would have no 
sense if one could not apply it to actually moving vortices. But in any case the requirement 
must be satisfied that characteristic velocity Uq of viscous vortex motion influenced by 
magnetic flelds comparable with low critical magnetic field Hd , must be significantly 
smaller than the speed of electromagnetic waves. 

The velocity scale uq can be naturally estimated as 

Uq = /le/X 

where 



is the self-energy per unit length of long straight-line vortex, e/Xis the characteristic scale 
of Lorentz force also related to unit length, and /j, is mobility of the vortex core. Below it 
will be seen that so defined Uq really serves as the velocity unit. If combine this definition 
and the known relations [4] 



with standard notations, an being the normal conductivity and A ~ 2kBTc being the 
order parameter, one obtains 

XhsTc Hci 

As a typical example, at Tc ~ 100 , A ~ 3 ■ 10^^ cm and Hc2/Hci ~ 500 , one gets 
the estimate Uq ~ 10^ cm/s . This value looks small enough to allow for applicability of 
quazi-static London approximation. In fact, such an approach was used in large number 
of works on motion of separate vortices as well as vortex lattices. The obvious exception 
is very dense lattice, with small inter-vortex distancies of order of coherence length. But 
our present subjects of interest are far from such complications. 
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III. EVOLUTION EQUATIONS 



In the framework of London approximation, the free energy E of vortex, placed 
into a given surroundings, is completely determined by the shape of its core, R{p) ~ 
{X{p),Y{p), Z(p)} , with X,Y,Z being coordinates of the core points and p being a 
scalar parameter. In accordance with the principles of mechanics and nonequilibrium 
thermodynamics, the simplest equation of a massless viscous evolution of the core line 
looks as 

pi-'dR/dt = f{R) (1) 

with Lorentz force on the right-hand side and friction force on the left, both being related 
to unit core length. 

By its sense, the parameter //"^ is the effective drag coefficient which is determined 
by all the dissipative energy losses conjugated with the core motion. Generally, there are 
at least two sorts of dissipative processes accompanying the motion (see, for example, 
the review [4]), namely, relaxation of the order parameter and normal currents induced 
by time-dependent own magnetic field of the vortex. A concrete expression for /i can be 
derived from more detailed theory, for instance, from the Ginzburg-Landau functional 
approach, under its reduction to London approximation [4]. The reduction is possible 
because normal self-current of moving vortex and corresponding dissipation are located 
mainly in a close vicinity of the core line, at distance comparable with coherence length. 
After the transition to London's description, the effect of normal currents becomes hidden 
in fj, , but these currents give no contribution to the Lorentz force [4]. Therefore, the 
reduction results in the identity whose meaning is balance of friction force and Lorentz 
force, as it is stated by the Eq.l, with f{R) being determined only by supercurrents. 

To write f{R) , one has not to evaluate the supercurrent distribution. Instead, as in 
general in mechanics and statistical thermodynamics, f{R) can be expressed by means 
of £"s variation under a small displacement of a local core fragment, that is as the func- 
tional derivative SE/SR{p) . However, the latter itself is not invariant with respect to 
arbitrary (non-degenerated) transformations of the parametrization R{p) and to physical 
dimensionality of p . In case of isotropic media, the only true invariant expression for the 
Lorentz force is 



^ dp 5E 



dL SR{p) 



OR 



dp 



6E $0 

- X AT] (2) 



SR{p) c 



Here the vector dR/ dp = R' is locally parallel to the core, dL = \R'\ dp is the differential of 
the core length, N = R' \R'\^^ = dR/ dL , and J is the density of full effective supercurrent 
which streams around core and pushes a given core fragment. 

The energy E = E{R{p)} includes self-interaction of vortex and its interaction with 
surroundings, in particular, with other vortices. Correspondingly, in general J consists of 
external currents and self- current of vortex determined by its distorsion. The Eqs.l and 
2 could be directly extended to a number of interacting vortices. Besides, in principle, 
one may add into E also interactions with pinning potentials. However, below we are 
interested only in motion of separate vortex in absence of pinning. 
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The parameter p enumerates strictly the core points. But in practice it is preferable 
to use another kind of parametrization, concretely, to introduce the parameter q which 
enumerates some suitable continuum of surfaces Q(r) = g , r = {x, y, z} , each possessing 
only one intersection with core line. The connection between p and new parameter q is 
implied by the obvious relation Q{R(p{q,t),t)) — q , and simple algebraic manipulations 
lead to the modified form of the evolution equations, 



dR 
'di 



OR dQ{R) 



dq 



dR 



dR 



dq 



SE 



(3) 



where the symbol (8> denotes the tensor product of two vectors. 

These equations describe how the intersection points marked by q move along the 
corresponding surfaces Q{r) = q ■ Clearly, this is factually two-dimensional motion. This 
feature becomes quite obvious if it is possible to identify q as one of cartezian coordinates, 
that is to use parallel planes as the marking surfaces. For instance, if thats are XY-planes, 
q — Z and Q{r) — z , then the Eqs.3 reduces to the equation 





/ 6E/6X{Z) \ 
[ SE/SY{Z) ) 



(4) 



with shortened notations X' = dX/dZ , Y' = dY/dZ . The Eq.4 describes the time 
evolution of X and Y coordinates of the core points marked with their Z-coordinate. 



IV. BOUNDARY CONDITIONS AND LOCAL APPROXIMATION 

In absence of pinning and more vortices, the vortex energy E = Eg + Ei consists of 
two parts: the energy Ei of the vortex interaction with transport or Meissner supercurent 
and the self-energy Eg. Therefore, the current in the Eq.2 also can be devided into two 
parts, J = Jg + Ji . 

Formally, Eg is a comphcated spatially non-local functional [5] depending on both the 
core configuration R{p) and the shape of sample. Among other factors. Eg includes the 
vortex interaction with the sample boundary what can be interpreted as attraction of the 
end fragments of the core to their mirror images placed outside superconductor. 

But, if the curvature radius of the core everywhere is not too small as compared 
with A , and besides, if the core nowhere is too close to itself, then the so-called local 
approximation is possible, 

Eg^sL = e j \dR{p)\ (5) 

where L is the core length. This well known approximation was argued and used as long 
ago as in 1968 by Galaiko [6], and later by many other authors (in particular, in [2-5]). 
Our own computer simulations showed that the relative error of evaluation of self-action 
force by means of local approximation does not exceed a few percents even if the curvature 
radius is as small as 0.1 A . 

In the local approximation the Eqs.l and 2 take the form 
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dR $n 
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dR 



dp 



(6) 



Here Jg is the self-current what flows through the very core. Its absolute value is inversely 
proportional to the local curvature radius of core. 

However, the local approximation needs to be accompanied by correct boundary con- 
ditions. The true conditions should take into account the vortex interaction with su- 
perconductor boundary. There are two ways to show that this interaction results in the 
orthogonality of the end fragments of core to the boundary. Thought these conditions are 
known at least since [6], sometimes thats are neglected, so it is desirable to present more 
argumentation. 

First, let us note that the force vector f{R) is always perpendicular to the local core 
direction. Indeed, any variarion 6R parallel to this core direction, 6R \\ dR/dp , merely is 
identical to a change of parametrization, without factual change of the shape, so it has no 
physical meaning and should result in 5E — (therefore the last expression in (2) always 
is consistent with previous ones). The same is seen from (6). As a consequense, any core 
point displaces perpendicularly to the core, in particular, the end points do which are 
placed just on the boundary. Hence, we must conclude that the end fragments always are 
oriented to be orthogonal with respect to the boundary. 

Secondly, the non-orthogonality would mean that the contour formed by core and 
its mirror image is broken at the end point, i.e. has infinitely small curvature radius 
here. From the point of view of exact Eg [5], if such a sharp "knee" occured it would 
cause infinitely strong Lorentz self-action force and consequently would be immediately 
straightened thus restoring the orthogonality. 

But, we must to underline that the principal conclusions to be deduced do not refer 
to the local approximation and can be derived from general non-local Eqs.l and 2 only. 



V. STICK-ELASTIC VORTEX TRANSFORMATION 
IN CURRENT-CARRYING PLATE 

To avoid a complicated mathematics, we confine ourselves by simplific superconductor 
geometry. Consider the vortex evolution in an infinitely wide plate, —D<Z<D, 
without pinning but in presence of transport surface supercurrent uniformly distributed 
over the boundary planes and obeying the London equation. If this current flows along 
Y-axis then 

j. = j. = o , j, = ^^H^,j{z) , m = h^^^^^ 

with h being the dimensionlcss measure of current density. 

Let initially the vortex pierces the plate in Z-direction being described with R = 
{X{Z, 0) = 0, 0, Z} . It has similar orientation soon after nucleation near the edge 
of a real finite plate. Then, due to obvious spatial symmetry, the vortex will remain 
inside the XZ-plane y = and keep only one intersection with any of XY-planes. In this 
situation the Eq.4 can be applied and, besides, reduced to only equation for X-coordinate, 
X — X{Z, , as a function of time and Z -coordinate: 
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The energy can be expressed as 



E — Eg + Ei — 



E. 



'^s-jlx{Z,t)j{Z)dZ 



(8) 



where the integral represents the energy Ei of vortex interaction with transport current 
(this expression differs only by some constant from the general E^ representation [5]). In 
the local approximation (5), the Eq.7 looks as strongly nonlinear diffusion- type equation 



with notation X" = d'^X/dZ'^ and characteristic velocity Mq introduced in Sec. 2. 

Here the left side is responsible for the friction, and two terms on the right-hand side 
represent the self-action force and transport current-induced Lorenz force, respectively. 
Clearly, because of the latter force both the vortex ends will forwardly move in one and 
the same direction (to opposite edge of the plate), while the middle of vortex will be more 
or less backward, and the larger is transport current the longer should be the distance 
AX = X{±D, t) — X{0, t) (below termed vortex stretching). 

Some predictions of further vortex behaviour can be deduced merely from the energy 
expression (8). Just after start the middle is still in rest. As the Eq.8 shows, in thick 
plate {D >> X) the unit displacement of every end leads to the Ei 's decrease by he . At 
the same time, the corresponding lengthening of each of two symmetrical core branches 
results in the Es 's increase by e per unit length. Consequently, if h > 1 then the total 
energy decreases, and the vortex stretchening along the drift direction becomes profitable. 
The lengthening process should last until the curvature of the most backward central part 
of the core becomes so large that the self-action force makes this part moving as quickly 
as the ends do. 

Thus, at h > 1 the vortex gets over the friction like clastic in water. Oppositely, 
at /i < 1 the stretchening is energetically unprofitable, and the vortex should move as 
deformed flexible stick. The transition from this stick-like behaviour to elastic-like one, 
when transport current increases from /i < 1 to /i > 1 , is the example of so-called 
" nonequilibrium phase transitions". 

Let us consider the steady drift of vortex as a whole, without change of shape. The 
corresponding solution on the Eqs. 7 or 9 can be written as X{Z,t) = ut + X{Z) . The 
stationary shape X{Z) and the drift velocity u — u{h, D) should be obtained from (7) or 
(9) with the help of above discussed orthogonality boundary conditions 




X 



+ ViTx^3{z) 



(9) 



dX 



{±D) = 



Besides, due to the mirror symmetry, the condition X'(±0) = should be satisfied. 

In this steady nonequihbrium state the self-energy Eg is constant, therefore, the work 
Mj produced by transport current per unit time. 



^ f'' J (z)dZ = 2t^/i£tanh(D/A) 

Ott C >/-D 

coinsides with the energy dissipation per unit time . In accordance with (1) and (2), 



Wa = - 



dR 



dt 



dL = ^J j QdZ = 2Su'^/fx 



where the notations 







are introduced. Wc took into account that actual displacement of the core always is locally 
perpendicular to its orientation. Only such displacements are physically meaningful and 
really cause the friction. Therefore, the drift velocity and the local core velocity are 
connected by the relation 



dR 



dt 



u 

Q 



Evidently, the factor Q determines at one and the same time local orientation of the core 

and degree of its stretching. 

Hence, the equality = Mj yields 

u hX ^ , ^ , ^ . hX 
U = — ^ — tanhm/A ^ — 10 
uq 6 d 

In view of above reasonings, at /i < 1 the vortex stretchening is weak, therefore, X'"^ 
is comparable with unit, Q ~ 1, AX ~ D and S ^ D. Then the Eq.lO shows that in this 
stick-like regime U ~ hX/6 ~ hX/D << h , i.e. the drift velocity is inversely proportional 
to the plate thickness. This is quite natural, because the surface current-induced Lorentz 
force acts only on the ends, while the friction almost equally acts on any core fragment. 

In general, the parameter 6 serves as the effective plate half-thickness. Obviously, 
always S < D . In the stretched elastic-like regime in thick plate anywhere at D — 
\Z\ » X the inequalities \X'\ >> 1 and Q >> 1 take place. Hence, S << D and 
what is essential it becomes almost insensitive to thickness. As a consequence, both 
the drift velocity and mobility u/h strongly increase as compared with stick- like regime 
and both become independent on thickness (below we shall see that 5 ~ A and C/ ~ /i 
, i.e. U becomes approximately D/X times larger). According to the M^s expression, 
the matter is that thought the energy dissipation Q times increases due to the core 
lengthening this effect is overpowered by its decrease because of Q times decrease of 
the factual core velocity ^ . As the result, the vortex stretching leads to smaller friction 
and smaller entropy production, under fixed vortex velocity, and to larger velocity ander 
fixed transport current. The picture looks as if most part of core slides along itself, but 
this process does not mean a real motion of core and so does not cause a friction and 
dissipation. 
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VI. DRIFT OF THE VORTEX ENDS 



To be convinced in what was said, let us consider vortex shape in the stretched elastic- 
hke regime. Because at AX » D most part of the core inevitably has a small curvature, 
it can be considered with neglecting the self-action. Then any of the Eqs.7 and 9 reduces 
to 

U ^ Vl + X'^j{Z) (11) 
Here from the characteristic exponential asymptotics does follow. 



exp 



A 



(12) 



Here Z > , X{—Z) = X{Z) , and for definitness the position X = is prescribed to 
the end fragments. It is easy to verify that corresponding self-action contribution in the 
Eq.9 indeed is negligibly small as compared with the current-induced force. 

We can get a rough estimate of the stretching if put on Z = in (12) and take into 
account that U > hX/D. Then the Eq.l2 yields 

AX/X -exp - > --exp - 



h "^KXJ D ^VA 

Hence, AX/X possesses exponentially strong dependence on D/X , and it can be giantly 
large if D exceeds A by an order of value or more. 

In view of this circumstance, the ratio AX/W with W being the width of a real finite 
plate, becomes of principal importance. Clearly, if AX » W then the steady drift of the 
vortex as a whole is impossible: the ends of vortex will achieve the opposite edge before 
the displacement of its backward central part will be comparable with W (all the more, 
before the velocity of this part becomes equal to that of the ends) . 

Consider the drift of the ends in such a non- stationary situation. Because the vortex 
lengthening is profitable, this drift can do independently on the motion of deepened 
backward part, as if thickness was infinitely large {D/X ^ oo) . To estimate the drift 
velocity, let us multiply the Eq.7 or 9 (with ^ =^ u) hy and integrate over variable 
z — D — Z from zero to infinity, with the condition X'{z oo) = oo which evidently 
corresponds to infinitely far backward center. Then both the Eqs.7 and 9 result in 



To evaluate this integral, note that in accordance with the orthogonality boundary 
conditions the shape of the end fragments of the core is parabolic, for instance, at upper 
end X{Z) = X{D) — {Z — DY /2p . with p being the curvature radius at the end point. It 
follows from the Eq.9 that X/ p = h~U . In this parabolic region the integration divergers 
but becomes cut after transition to exponential asymptotics (12). The estimate of the 
integral leads to approximate equation 

U^h{h- 1) /{h-l + \n[2zo{h - U) /A]} (13) 
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where zq is the depth of the crossover point, zq ^ AX . The Eq.l2 helps to estimate the 
end drift velocity in thick plate. Obviously, it turns into zero ai h ^ 0, in agreement 
with D — > oo limit of the estimate for stick-like regime. It can be shown that velocity of 
the steady drift of the vortex as a whole is only slightly smaller differing by a multiplier 
of order of unit. 



VII. GIANT VORTEX STRETCHING 

The exponentially large vortex stretching is the most significant possibility of vortex 
evolution in thick plate, as well as in bulk samples in general. Therefore it would be useful 
to more correctly justify the above simplific estimate of AX . 

Note that AX > L/2 — D . Divide both sides of (9) by j{Z) and integrate from zero 
to D . This results in 



/dZ r 
— — , A= arctan(X') 



d A 



dZj{Z) 



dZ 



It is easy to notice that 



4 A 



m m 



so AX > B — Aq — D. The calculation of integral B gives 

B = ^cosh(L>/A){arctan[exp(L>/A)] - j} 
Then, after simplifications possible due to D >> A, one finally obtains 

AX/X > ""^^ ~ exp{D/\) 



(14) 



with A; < 1 . Because u is monotonously growing function of h , the coefficient in front 
of exponent is positive if h exceeds some level larger than unit, for example, if /i > 2 . 
Hence, ate least at h > 2 the vortex stretching is exponentially strong. 

This estimate is obtained in the framework of local approximation. More correct 
estimate should give a lesser value, because of self-attraction of the core in middle part 
of the plate where two symmetrical exponential tails described by (12) meet one another 
and form an arc. Such a non-local effect is most essential just under the specific plate 
geometry. However, the non-local correction can not change the shape of the front vortex 
part where the non-local interaction is weak as compared with other forces. It is not hard 
to show that the latter requirement is satisfied if \Z\ > Zq, where Zq is the solution on 
equation 



hex-p[-{D - \Z\)/X] ^ VA/27r \Z\ exp(-2 \Z\ /A) 

If take into account that the more is h the less is Zq, then this equation yields Zq < D/3 
. Hence, at \Z\ > D/3 the mutual attraction of two core branches can be neglected, and 
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the asymptotics (12) remains valid. This means that the maximaUy possible effect of 
non- locality is the replacing D in the exponent by qD with a > 2/3 . Consequently, the 
lower bound for the stretching with confidence can be estimated as 

AX/A > exp(2L)/3A) 

Thus, even in the worst case the non-local effects do not abolish the exponential character 
of stretching. 

For example, if A ~ 3 • 10~^ cm and h equals to a few units, then even at ~ 20A < 
10~^ cm one gets /S.X > mXexp{2D/3X) , with m ~ 1 , i.e. AX > 1cm what exceeds a 
width of any realistic sample. Thus at first the vortex should form a ring whose shape 
approximately copies that of the sample cross-section. During this process the velocity of 
backward deepened core part is primarity determined by its distorsion which is created 
in the beginning of stretchening and thus has curvature radius of order of D . Hence, this 
velocity is of order of UqX/D , and at the moment when the ends will meet one another 
the displacement of most backward point will be yet as small as ~ XW/Dh «W . 

VIII. DISCUSSION AND RESUME 

It seems clear that both the above conclusions can be extended to bulk current-carrying 
superconductors with another geometry, for instance, to round wires, if treat 2D and W 
as minimal and maximal diameters of cross-section of the wire, respectively. 

Due to possibility of giant deformation and stretching of vortices, the thermodynami- 
cally nonequilibrium process of vortex penetration can promote formation of complicated 
many-vortex dynamical configurations which seem rather strange and unprofitable from 
the point of view of equilibrium thermodynamics. The presence of an external magnetic 
field parallel to transport current should lead to formation of spiral-like configuration 
instead of ring-hke one and thus especially ensure the entangling of vortices. The Eq.4 
enables us to describe this scenario in details, if choose Z-axis to be directed along the 
wire. Besides, the presence of weak pinning should amplifier the stretching of vortex and 
additionly complicate its shaping, because the motion of deepened part of vortex is char- 
acterized by relatively small forces of order oi eX/D (much smaller than forces ~ e what 
act on the end fragments) and so may be easily held back by pinning centers. 

We would like to underline the role of orthogonality boundary conditions. In the work 
[3] the equation was under use similar to our Eq.9, but boundary conditions was formu- 
lated in terms of the tension of core line. One can see from [3] that such conditions make 
it impossible to consider the case of high surface transport current > ifci corresponding 
to the elastic-like regime. 

To resume, we formulated the invariant equations of viscous motion of arbitrarily 
shaped 3D vortex lines, and applied them to careful analysis of the scenario of vortex 
penetration into a thick superconducting sample. As it was argued, the vortex-vortex 
interaction does not significantly affect the penetration process. But, of course, a full 
description of resistive state leads to more complicated tasks about vortex-vortex inter- 
actions deep inside the sample. 
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